Growing condensate in two-dimensional turbulence 
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We report a numerical study, supplemented by phenomenological explanations, of "energy con- 
densation" in forced 2D turbulence in a biperiodic box. Condensation is a finite size effect which 
occurs after the standard inverse cascade reaches the size of the system. It leads to emergence of 
a coherent vortex dipole. We show that the time growth of the dipole is self-similar, and it con- 
tains most of the injected energy, thus resulting in an energy spectrum which is markedly steeper 
than the standard k~^^^ one. Once the coherent component is subtracted, however, the remaining 
fluctuations have a spectrum close to . The fluctuations decay slowly as the coherent part grows. 

PACS numbers: 47.27.E-,92.60.hk 



A big difference between 2D and 3D turbulence is the 
generation of large scale structures from small scale mo- 
tions pj, 121]. This occurs because, if pumped at interme- 
diate scales, the 2D Navier-Stokes equations favor en- 
ergy transfer to larger scales [1, 0, HI, 0], a phenomenon 
known as an inverse cascade. Simulations [?, "s] and ex- 



periments [9l, lid, UJj show that large scale accumulation 



of energy is observed if conditions permit the energy to 
reach the system size. In this letter, we study the "con- 
densate" emerging in the form of two coherent vortices 
in a biperiodic box in 2D. Let us begin by briefly re- 
viewing the classical 2D turbulence theory of Kraichnan, 
Leith and Batchelor (KLB) [i, 0, The essential dif- 
ference with 3D turbulence is the presence of a second 
inviscid invariant, in addition to energy, the enstrophy. 
Stirring the 2D flow leads to emergence of two cascades. 
Enstrophy cascades from the forcing scale, /, to smaller 
scales (direct cascade) while energy cascades from the 
forcing scale to larger scales (inverse cascade). Viscosity 
dissipates enstrophy at the Kolmogorov scale, r^, which is 
much smaller than / when the Reynolds number is large. 
The energy cascade is blocked at a scale C ^ ^ by 
a frictional dissipation (usually due to friction between 
the fluid and substrate although other mechanisms can 
be imagined) after a transient in time quasi-stationary 
regime. Then a stationary KLB turbulence is established 

3,13, Q- Applying Kolmogorov phenomenology (see e.g. 

l|]) KLB predicts an energy spectrum scaling as in 
the direct cascade, and as k~^^^ in the inverse cascade. 
Here k is the modulus of the wave- vector. The KLB 
spectra imply that velocity fluctuations at a scale r, 6vr 
scale as e^/^/~^/^r and (er)^/^ in the direct and inverse 
cascade ranges respectively. KLB theory is confirmed by 
simulations [12, 13] and experiments [21, fllj, where a suf- 
ficient range of scales was available to form the cascades. 
If the frictional dissipation is weak so that ( exceeds the 
system size L then ultimately the "condensate" regime 
emerges (sl, where the standard KLB does not apply. 

One of the primary motivations for studying 2D tur- 
bulence is that it is structurally and phenomenologically 



similar to quasi-geostrophic turbulence [l^] which de- 
scribes planetary atmospheres [l3]- In addition, a recent 
resurgence in theoretical interest in the condensate state 
was sparked by experimental [11] and numerical [18] ob- 
servations of large scale coherent vortices associated with 
energy condensation in forced, bounded flows. In this let- 
ter, we report results of a set of numerical experiments 
designed to give a clean, detailed study of this energy 
condensation phenomenon in its own right. 

We solved the incompressible forced Navier-Stokes 
equations with hyperviscous dissipation in 2D: 



dfU + {u ■ V)u + Vj9 = uA^u ■ 
V-n = 0. 



(1) 



The domain is a doubly-periodic box of size, L = 27r. 
The forcing, /, injects velocity fluctuations and energy 
at an intermediate scale / with energy injection rate e. 
For the simulations shown in Figs. [T] and [3l / = 27r/50 
and e = 0.004. We use a standard pseudo-spectral solver 
with full dealiasing. The resolution varied from 256^ to 
1024^. For developed condensate computations, the res- 
olution was only 256^ owing to the requirement of in- 
tegration for tens of thousands of forcing times, which 
was done using a 3rd order Runge-Kutta integrator with 
integrating factors. The timestep was decreased as the 
condensate grows such that it satisfies At < coAx/umax, 
where Ax is the grid spacing, iimax is the maximum ve- 
locity and Co is conservatively taken in the range 0.2 - 
0.5. Energy injection was done in a spectral band us- 
ing a stochastic additive force with fixed amplitude and 
random phase. The correlation time is the numerical 
timestep. Small scale dissipation was provided by 
hyperviscosity which is not expected to affect the large 
scale behavior. There is evidence [19], and our additional 
simulations (not discussed here) are in support, that the 
presence of large scale damping can considerably compli- 
cate the description of the large scales. 

The forcing described above is short correlated in time 
and characterized by the energy injection rate, e, and the 
forcing scale, /. The majority of the injected enstrophy 
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FIG. 1: Vorticity snapshots at a succession of times : (A) t = 100, (B) t = 2000, (C) t = 10000. Decomposition of the t = 10000 
snapshot into coherent, (D), and fluctuating parts, (E). Energy spectra of the full field, coherent part and fluctuating part, (F). 



cascades towards smaller scales to be dissipated by vis- 
cosity. In our numerics, the closeness of to / meant that 
only 40% of the injected energy goes upscale from / with 
the rest going downscale with the enstrophy. We simu- 
late the zero friction case to assure that eventually the 
energy starts piling up at the system size scale. The di- 
rect cascade is set up in a short time, ri ~ /^/^/e^/^. The 
time, r*, for the inverse cascade to populate all scales up 
to L, is much longer. Based on Kolmogorov arguments, 
r* ~ e~^/^L^/^. In the simulations, ^ 1000, in units 
where ri is about 1. 

At t > we observed a condensate consisting of two 
big vortices of the size ~ L separated by a hyperbolic 
domain of comparable size. Fig. [T] (A),(B) and (C) il- 
lustrates the phenomenon with a series of vorticity snap- 
shots. The condensate is formed to ensure that (a) the 
integral vorticity is zero in accordance with zero integral 
vorticity injected by the small-scale pumping, and (b) 
majority of energy brought by the inverse cascade is ac- 
cumulated at the largest scale, L. Two identical vortices 
rotating in opposite directions satisfy these conditions. 
Due to biperiodicity. Fig. [T] actually depicts the emer- 
gence of a vortex crystal. Such crystals have been ob- 
served both numerically [8] and experimentally [10]. The 
vortices drift slowly over time but the square symmetry 
of the crystal is preserved by this drift. 

Evolution of the vortices is slow relative to the back- 
ground fluctuations which permits a separation of the 
flow into coherent and fluctuating components in the 
spirit of 23,[21|. The highest amplitude coefficients of the 
wavelet transformed vorticity are assigned to the coher- 
ent component and the remainder to the incoherent com- 
ponent. Inverse wavelet transforms are then taken. This 



decomposition is shown in Fig. [T] (D) and (E). Fig. [TjE) 
shows that the fluctuating part is almost statistically ho- 
mogeneous, whereas the coherent part is strongly inho- 
mogeneous. Note that the characteristic amplitude of 
the vorticity fluctuations is larger than the coherent part 
of the vorticity over most of the domain. Ultimately we 
expect the coherent flow to dominate the fluctuations ev- 
erywhere but we have not reached this regime. 

As seen in Fig. [2](A), one observes oc growth of the 
maximum value of the coherent part of vorticity with 
time. Furthermore, simulations show the global growth 
oc of the coherent velocity profile. This global self- 
similarity is evident from Figs. [2](B) andlSfC). The law 
(X is naturally explained by the energy accumulation 
injected at the constant rate, e, by forcing. In the hyper- 
bolic region one estimates the coherent velocity as y/ei. 

The mean velocity profile within the vortex is almost 
perfectly circular. To a good precision higher order har- 
monics are suppressed relative to the zeroth order one. 
The velocity profile deduced from the simulations fits, 
is (X r~^, where ^ ^ 0.25, in the range, L ^ r ^ 
and thus the vortex core is roughly /. This is illus- 
trated in Fig. [21(D) showing the equivalent vorticity pro- 
file, (X We plot r~^-'^^ profiles for two different 
forcing scales to check that the profile is insensitive to it. 

So far in this letter we have been discussing the spatio- 
temporal features of the condensate. Complementarily, 
and following tradition developed in turbulent studies, 
one may also analyze velocity spectra. Evolution of the 
spectra in time is shown in Fig. (|3]), where one clearly 
see transition at from k~^^^ to scaling steeper than 
that is numerically close to . Similar state- 
ments were made before in Refs. 22, ^3|. We claim 
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FIG. 2: Self-similar growth of the condensate. (A) Maximum 
vorticity as a function of time, (B) Angle-averaged vortic- 
ity, r2(r), as a function of distance, r, from the vortex center 
for successive times, (C) Same profiles rescaled by a/^ — t*, 
(D) Q(r) in the developed condensate regime for 256x256 and 
1024x1024 simulations with two different forcing scales. 




FIG. 3: Time evolution of the spectrum showing the transi- 
tion from a standard k~^^^ inverse cascade to a scaling 
at large scales. Spectra have been locally averaged in time. 



that this spectrum is spurious and should not be taken 
as evidence of a cascade in the KLB sense. The coherent 
part of the flow has almost no fluctuations and, if it is 
removed, the steeper than k~^^^ scaling disappears en- 
tirely. By contrast, the enstrophy cascade of KLB 
involves fluctuating vortices across many scales. 

The spectrum of the fluctuations is shown in Fig.fljF), 
it is close to k~^. There is some evidence that shortly 
after r^^ the fluctuations retain a k~^^^ spectrum as was 
observed in [22|,[23|]. However, this regime is transient and 
does not persist for long after the condensate emerges. 



Larger resolution will be required to settle this ques- 
tion unambiguously. Fluctuations play a relatively mi- 
nor role in the overall energy balance with the majority 
of the energy absorbed into the condensate as shown in 
Fig. (|4])(A). They contain more of the enstrophy as shown 
in Fig. (|4])(B), but they decay in amplitude as the con- 
densate grows so that the flow becomes more and more 
coherent as time passes. The data suggest that decay of 
the amplitude of the background fluctuations is logarith- 
mic, or very weakly power law. 

We now present an attempt to describe phenomeno- 
logically, the universal nature of the asymptotic conden- 
sate state. Consider an individual vortex at t ^ r*. 
It has a core of radius ~ / and its spatial extent is es- 
timated by the system size, L. In discussing the spa- 
tial structure of the vortex, for example its mean vor- 
ticity profile, Q = (V x n), we will track its depen- 
dence on the distance, r, from the center of the vortex, 
Q{r). Once the almost circular vortex emerges, it keeps 
sucking energy from the turbulent background which can 
be approximately described by an inhomogeneous eddy 
diffusivity, D{r). Another large scale characteristic af- 
fected by the eddy-diffusivity is the fluctuation enstro- 
phy, H{r) = ((V X n — 1])^). The focused regime is adia- 
batic so that the equations governing the quasi-stationary 
radially symmetric distribution of ft and H on the top of 
the turbulent background are the eddy-diffusivity equa- 
tions 



drvDdrVt = 0, drvDdrH = 0. 



(2) 



Naturally, D can be expressed in terms of the typi- 



cal Lyapunov exponent. A, D 



"A. In homogene- 



neous turbulence A would be self-consistently estimated 
as ~ V^. However present situation is inhomogeneous, 
with a strong, ~ shear. Mixing in the presence of 
strong shear was discussed in [24] . It was shown that the 
dependence of the effective Lyapunov exponent on the 
mean shear, ~ 1], and the background enstrophy, can 
be estimated as 



A 



(3) 



Here th is the correlation time of the background vortic- 
ity fluctuations. The actual nonpar ametric regime we 
are interested in is thX ^ 1. Thus keeping the two 
asymptotics in Eq. (|3]) will, in fact, give upper and lower 
bounds. Returning to Eqs. ([2]) one notes that the phys- 
ically meaningful solution for Q corresponds to a flux 
state zero mode of the eddy-diffusivity equation, describ- 
ing accumulation of the mean vorticity (and the mean 
energy) inside the vortex, Drdr^ = const (wrt r). On 
the contrary, the physically meaningful solution of the 
eddy-diffusivity equation for H is the one corresponding 
to a zero spatial flux, H = const (wrt r). Note, that this 
spatially homogeneous distribution of H is in agreement 
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with the results of simulations. Combining all these esti- 
mations with the global energy conservation one arrives 
at the following bounds 



ri/A<Cl; 
thX > 1. 



(4) 



These estimates for the mean vorticity profile are fully 
consistent, in describing both the overall temporal dy- 
namics and the exponent of the mean vorticity profile, 
with the aforementioned numerical simulations, shown 
in Fig. 2: 1/5 < ^ 0.25 < 1/3. The corresponding es- 
timate for the spatially homogeneous enstrophy, also ex- 
pressing direct enstrophy balance at the pumping scale, 
is H ^ £/{J?\). This formula, combined with Eq.([3]) for 
the Lyapunov exponent, predicts a slow algebraic decay 
of the background enstrophy in time. This is again con- 
sistent with simulations. 

Finally, the spatial homogeneity of H suggests that 
majority of the injected enstrophy cascades to smaller 
scales, r <^l. However, a subdominant portion will also 
penetrate to the larger scales, e.g. resulting in the 
spectrum observed in the simulations, see Fig. 3. An 
explanation of the spectrum observed at the scales 
larger than / after subtraction of the coherent component 
of the flow is as follows. In a range of scales smaller than 
L vorticity fluctuations are advected passively. Passive 
scalar theory, developed in [iH, predicts an ~ 1/r^ de- 
cay for the pair correlation of a scalar at the scales larger 
than injection scale in two dimensions and for non-zero 
value of the Corrsin invariant, which is the integral of the 
pair correlation function of the pumping. However, vor- 
ticity is a curl of velocity injection and thus the vorticity 
is injected at / with zero value of the Corrsin invariant. 
This leads to the localized, ~ expression for the pair 
correlation function of vorticity, which in turns translates 
into the observed spectrum. Notice that similar ex- 
planation for this scaling, referred to as passive inverse 
energy cascade, was reported in [2^. 

To conclude, we performed numerical simulations of 
energy condensation in forced 2D turbulence. We split 
the flow into coherent and fluctuating parts, observed the 
power-like shape of the coherent vortices and self-similar 
growth in time of the coherent flow. The coherent struc- 
ture is responsible for the spurious energy spectrum 
observed in previous numerical experiments. The fluctu- 
ations have an energy spectrum oik~^ and they diminish 
in importance as the condensate grows. We also pre- 
sented phenomenological description of the simulations. 
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clear Security Administration of the U.S. Department of 
Energy at Los Alamos National Laboratory under Con- 
tract No. DE-AC52-06NA25396. 
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FIG. 4: (A) Time evolution of the energy contained in the 
condensate and background fluctuations (B) enstrophy. 
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